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Abstract 

Gauge fixing in the non-perturbative domain of non-Abelian gauge theories 
is obstructed by the Gribov-Singer ambiguity. To compare results from dif- 
ferent methods it is necessary to resolve this ambiguity explicitly. Such a 
resolution is proposed using conditions on correlation functions for a fam- 
ily of non-perturbative Landau gauges. As a consequence, the various re- 
sults available for correlation functions could possibly correspond to differ- 
ent non-perturbative Landau gauges, discriminated by an additional non- 
perturbative gauge parameter. The proposal, the necessary assumptions, 
and evidence from lattice gauge theory calculations, are presented in detail. 
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1. Introduction 

To describe gauge-dependent degrees of freedom, like quarks and gluons, 
it is necessary to fix a gauge. In the non-perturbative domain of non-Abelian 
gauge theories this is complicated. Local gauge conditions, e. g. for Landau 
gauge 

PAl = 0, (1) 

are satisfied by more than one field configuration A^, being Gribov copies of 
each other Jl|. Conditions, which are not local, are required to resolve this 
ambiguity [2j. 

It is desirable that such supplemental conditions are formulated indepen- 
dent of the particular algorithmic implementation of a method. Furthermore, 
in the context of lattice calculations and functional methods, a formulation 



using correlation functions would be desirable. In the perturbative case the 
latter is possible. E. g., ([I]) can be formulated as the vanishing of the longi- 
tudinal part of the gluon propagator D^, 



Also, certain supplemental conditions for the Landau gauge can be formu- 
lated in this form. One example is the absolute Landau gauge, which requires 
to absolutely minimize the total trace of the gluon propagator 0, 0] 



with c a positive constant. However, ([3]) requires regularization in the con- 
tinuum and it is not clear whether this introduces new ambiguities. Thus, 
alternatives are desirable. This should in general be possible by conditions 
on correlation functions: Given a set (of configurations) of sets of distinct 
events (Gribov copies identified by the gauge field at every space-time point), 
it is always possible to construct probability distributions (gauges selecting 
a particular Gribov copy) such that they are distinguished and identified 
uniquely by at least some of their - finite or infinite - moments (correlation 
functions). The remaining question is whether this possibility can be cast 
into a practical method, in particular whether the distinguishing correlation 
functions can be determined. 

In continuum studies a one-parameter family of correlation functions has 
been obtained 0, 0], distinguished by the ghost propagator, an infinite- 
order moment with respect to the elementary gluon field. This motivates 
the present study whether this family could indeed arise as the consequence 
of a gauge choice, as speculated in 0,0]. Indeed, at least in a finite volume a 
corresponding family is constructed using gauge conditions resolving, at least 
partly, the Gribov-Singe ambiguity. Assuming the (qualitative) correctness of 
the results in the continuum this would suggest a one-to-one correspondence. 

2. The structure of the residual gauge orbit 

After imposing the Landau gauge condition (JTJ remains a set, called the 
residual gauge orbit j3[ here, of Gribov copies separated by large gauge trans- 
formations. These orbits can be divided in Gribov regions [if, with the orbit 
of each configuration passing at least once through each region [7]. A pos- 
sible first step to construct a definite gauge-fixing prescription is to restrict 



Pll p u D^ = 0. 



(2) 




(3) 
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Figure 1: The average number of Gribov copies as a function of lattice extension (left panel) 
and volume (right panel) in two, three, and four dimensions at a — 0.22 fm, distinguished 
by their value of trD. At finer a at fixed physical volume the number of copies increases 
further [lij . All lines drawn to guide the eye. Lattice volumes N d are for d = 2, 3, and 
4 from the sets {10,18,26,34}, {8,14,20,26}, and {6,10,14}, respectively, throughout. 
A fit of type Aexp(V/e d ) gives an e of approximately 16, 4.4, and 2.9 fm in 2, 3, and 4 
dimensions, respectively. 



the residual orbit to one of these regions. Here, the first Gribov region where 
the Euclidean Faddeev-Popov operator —d^D^, with the covariant derivative 
Dp, is positive semi-definite is chosen. 

For this restriction an explicit prescription purely in terms of Green's 
functions is not yet known. However, general considerations for a positive 
operator and all results so far (for a recent compilation of references see 0, 0]) 
strongly suggest that it is sufficient that the ghost propagator D G , being the 
expectation value of the inverse Faddeev-Popov operator, is required to be of 
definite (negative) sign and possibly monotonous. An explicit evaluation in 
1+1 -dimensional Coulomb gauge [8] is finding exactly that the ghost propa- 
gator changes sign for some momenta when evaluated outside the first Gri- 
bov region. However, lattice calculations by construction jsj and functional 
methods implicitly lOj] ensure that results are obtained from inside the first 
region, and this question is of minor concern. Hence, it is not necessary to 
obtain the spectrum of the Faddeev-Popov operator explicitly. 

Still, inside the first Gribov region there are many Gribov copies [4]. On 
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Figure 2: The distribution of trD compared to 6(280 MeV, oo) on a 20 3 lattice at a = 0.22 
fm. The left panel shows for an example configuration that for each distinct Gribov copy, 
identified by differing trD values, also the b values differ. 22 copies have been generated 
for this configuration, with fO distinct Gribov copies found, one peak being hidden by 
foreground peaks. The right panel shows the distribution for f450 configurations. 



a finite lattice^, it appears that Gribov copies can be distinguished by their 
difference in the value of the trace of the gluon propagator ([3]) jl, [ll| . This 
will be assumed henceforth. The number of Gribov copies rises quickly, see 
figure [H and thus there is a dense, though not necessarily connected, set of 
them in the infinite- volume limit@. Here, it is investigated whether other dis- 
tinctions between Gribov copies exist. Motivated by the continuum studies 
0, 0] the ghost propagator was evaluated. It has been checked whether the 
Gribov copies can be characterized by the renormalization-group invariant 
quant it^ 

to, p) = m = < ?D ^ (4 ) 

{P - ' G(P) P 2 D a (P) W 

for fixed p ^ P and G(p) = —p 2 Dg{p) is the ghost dressing function. If 
both momenta in the definition (jlj) are large and comparable, b will be of 



1 For details of the lattice calculations see 0, [t| and for the generation of different 
Gribov copies fill ] . 

2 In two dimensions for the volumes investigated in |12j only a negligible amount of Gri- 
bov copies are present 0]. In these volumes a so-called scaling- type behavior 0] prevails. 

3 This definition is prone to violations of rotational symmetry. For simplicity here all 
ghost momenta are evaluated along the x-axis, but a more robust method is required 13]. 
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order one up to sub-leading corrections. From a practical point of view it 
is therefore better to use p small and P large. Here, p will be taken to 
be the smallest accessible momentum, zero in the infinite- volume limit, and 
P the renormalization scale p taken in the perturbative domain, which can 
be chosen infinite in dimensions smaller than 4. In that case, G(P) = 1 
trivially, and the normalization can be dropped altogether. Note that the 
choice directly influences the severity of finite-volume artifacts, as they are 



stronger the smaller p is chosen. Other choices will be studied elsewhere [13 

Indeed, Gribov copies with differing trD are found to have also differing 
values of b, as can be seen for an example configuration in the left panel 
of figure [2J The quantity b measures in a sense the proximity to the Gribov 
horizon to the extent that it is dominated by the lowest eigenvalue and eigen- 
state of the Faddeev-Popov operator. It is found that tr.D and b turn out 
to be almost uncorrelated, see figure |2] right panel. This would imply that 
there is little relation between the values of b and the so-called fundamental 
modular region of minimal trD values. However, it cannot be excluded that 
this is a lattice (discretization or finite volume) artifact. 

2.1. Constructing a gauge 

Since b appears to be different for each Gribov copy, at least on the vol- 
umes investigated here, this would be sufficient to construct a non-perturbative 
extension of the Landau gauge, called Landau- B gauges hereafter. It is de- 
fined by requiring that the function 6(0, p) (6(minp, p) on the lattice) should 
have a prescribed value of B on the average. This is implemented by the 
following prescription: Take from each residual gauge orbit the copy with b 
closest to B, and average over these copies to obtain the gauge-fixed quanti- 
ties. In case of exactly the same difference for two copies, select one randomly 
to ensure the correct value of B on the average. By construction, the ghost 
propagator will then satisfy b = B on the average, if the value of B lies within 
the range of possible values for b. Otherwise a value as close as possible to 
B will be obtained. 

This prescription is valid for any volume, and specializes the perturbative 
Landau gauge ([[]) if more than one possible value of b exists. If each Gribov 
copy of a given orbit has a unique value of b, this actually resolves the Gribov 
ambiguity completely. If not, this implies that further additional conditions 
can be imposed to lift the remaining degeneracies, specializing the Landau- 
B gauges defined here further. If, like here, this is not done and there are 
degeneracies, just a random representative among the remaining degenerate 
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Gribov copies is chosen. This is the same as done in minimal Landau gauge 
(sf, where not even different values of b are distinguished. 

Note that both, Landau- B gauges and minimal Landau gauge, are not af- 
fected by any Gribov ambiguity anymore: Without specifying how to resolve 
potential further degeneracies they yield a result which gives an expectation 
value over degenerate Gribov copies. In a thermodynamic setting this would 
be the equilibrium result, i. e. the most likely one, for correlation functions, 
provided the implemented algorithm faithfully represents the distribution of 
Gribov copies. In particular, no Gribov problem, i. e. the necessity to know 
all Gribov copies, arises in general. This problem is present, e. g., in the 
absolute Landau gauge or if the most extreme values of B are desired. 

Here 'equilibrium' is not indicatin g th e existence of a preferred value: 
In the sense of stochastic quantization [lOj , no preferred Gribov copy on the 
orbit exists, any choice is equally admissible. Specifying a unique copy would 
therefore imply a maximum non-equilibrium situation, which is nonetheless 
as valid a choice as the equilibrium value of the minimal Landau gauge. 

2.2. Connecting to functional equations 

If this indeed identifies a gauge also in the infinite-volume and contin- 
uum limit this implies that the value of the ghost propagator distinguishes 
different treatments of Gribov copies. By analogy imposing conditions on 
the ghost propagator as boundary conditions on the functional equations, 
as detailed in |6j, would then be equivalent to treat Gribov copies in such 
calculations. This would give an interpretation for the fact that the value of 
B cannot be determined self-consistently inside functional equations [6]: B 
has the meaning of a second non-perturbative gauge parameter, besides the 
conventional gauge parameter £ to distinguish covariant gauges. This pre- 
scription is not unique, for alternatives see jil Note that choosing B to 
be the average value among Gribov copies then provides a prescription how 
to implement the minimal Landau gauge in functional methods. If further 
degeneracies would exist, this interpretation would imply that with improved 
truncations compared to j(| further undetermined boundary conditions must 
appear in the functional equations to permit resolving these degeneracies. 
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Figure 3: The gauge corridor for different lattice extensions (left panel) and volumes (right 
panel). For each volume < min6(minp, fi) > and < max6(minp, (i) > are shown, where 
minp is the smallest non-zero momentum on each lattice, and fi is set to oo for two and 
three dimensions, yielding G(p) — 1, and to fi = 2 GeV in four dimensions. 



3. Results for the correlation functions 

The value of B can only be chosen in a gauge corridor of values which 
are assumed on the residual gauge orbit. This corridoi@ is shown in figure 
EJ Of course, if the corridor should decay into several bands with forbidden 
ranges in between, which is not observed here, choosing B in between bands 
yields possibly ambiguities. In three and four dimensions the lower bound of 
the corridor is almost independent of volume, and by construction bounded 
from below. The upper bound is strongly increasing with volume for those 
investigated here. In two dimensions, due to the near-lack of Gribov copies 
at these volumes, the effect is rather small. Whether the up per boundary 



diverges in the infinite- volume limit remains to be investigated [13| . Without 
explicit exploration it cannot be excluded that the corridor closes again at 
much larger volumes, corresponding to total degeneracies of Gribov copies 
with respect to b. 

To compare extreme cases, the correlation functions for the cases with 
average, and minimum and maximum attainable value of B, corresponding 



4 Inherently, this method to identify Gribov copies llfl cannot guarantee to find all 
copies. Hence, the boundaries of the corridor represent upper bounds for the lower bound- 
ary and lower bounds for the upper boundary. Finding the exact ones is a Gribov problem. 
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Figure 4: The gluon propagator (left panel) and the ghost dressing function (right panel) 
for various gauges from a 20 3 lattice at a = 0.22 ftn. The minimal Landau gauge is 
described in and the absolute Landau gauge in The minS and maxB gauges 
select the minimum and maximum possible value of B, respectively, and are described in 
the text. The inverse Landau gauge maximizes ([3]) over the set of all local minima on the 
residual gauge orbit [HI]. 



to setting B to zero and infinity, respectively, are determined. The results for 
the gluon propagator and the ghost dressing function in these three types of 
Landau- B gauges, compared also to gauges based on ([3l), are shown in figure 
HI The ghost propagator differs very strongly at momenta below 1 GeV for 
the three Landau- B gauges. In fact, the maxi?-case is enormously enhanced 
compared to the result m minimal Landau gauge! The gluon propagator 
is, at least for this lattice size, almost indistinguishable for all gauges. This 
could have been anticipated due to the lack of correlation of the ghost and 
gluon propagator gauge parameters, shown in figure El Arguments exist 13 



that it, and any finite moments, may coincide for any fixing of the residual 
gauge degree of freedom in the infinite volume limit. 

As the experience shows (see e. g. [3]), rather large volumes will be 
necessary to be able to extrapolate to infinite volume. This will be addressed 



in the future [131 ] . However, there are in principle three qualitatively different 
outcomes, if no further degeneracies exist. First, the corridor is of finite 
width with finite boundaries. Then a family of decoupling-type correlation 



5 Note that the bounds of [14( are trivially fulfilled for any £>-gaugc, as they depend on 
the chosen B-gauge. 
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functions is obtained. Second, the upper boundary is infinite, and the ghost 
propagator corresponding to 1/B = behaves as in the scaling case. The 
behavior of the gluon propagator is then not yet fixed, see the next subsection. 
Finally, the possibility exists that the upper bound is infinite, but the ghost 
propagator is not that of the scaling solution. If additional degeneracies 
exist, other possibilities are possible. E. g., the 6-corridor can close to a 
single value, finite or infinite. In this case, only one .B-gauge survives, with 
a unique ghost and gluon propagator. 

3.1. Relation to functional results 

Assume for the moment the correctness of the functional results yielding 
a one-parameter family of solutions with the scaling solution as an endpoint 
0| . Also an infrared finite gluon propagator is not at odds with this, as 



functional studies imply such a behavior at small volumes 16j . Even if the 
gluon propagator would remain infrared finite in the infinite-volume limit 
(and possibly S-independent), this would not be at odds with this, since this 



(corresponding to a critical exponent 17| of 1/2) is still compatible with a 



scaling-type behavior in case of mild angular variations of the ghost-gluon 
vertex [lT]]. This is compatible with existing results 0, 12, 18], and is expected 
from functional studies [l9j]. The role of the absolute Landau gauge is then 
not clear, as it so far appears unrelated to the 5-gauges. In particular, it 
could well be that the interpretation given in [3[ is wrong, and the absolute 
Landau gauge is not connected to a scaling behavior. This requires more 



investigations [13 . 



4. Summary 

Thus, using lattice gauge theory an one-parameter family of correlation 
functions at finite volume is found. These are distinguished by a second 
gauge parameter B imposed on the ghost propagator. If there is no further 
freedom in choosing a Gribov copyj, this would resolve the Gribov-Singer 
ambiguity completely. In addition, this would be achieved by conditions on 
the correlation functions only. 

As a consequence, even the qualitative infrared behavior of correlation 
functions could be determined by a gauge choice, and any physics contained 



6 In the case that further degeneracies exist, there would be even more freedom to vary 
correlation functions by choice of Gribov copies. 



9 



in there would be mixed with gauge contributions. A fuller analysis of this 



question will require more detailed studies [13], and in particular much larger 
volumes. This applies also to the connection to a corresponding family of 
correlation functions obtained in functional studies [6| , which could here only 
be established at the investigated volumes. If this could also be established 
for larger volumes and eventually in the infinite- volume limit, this would 
be extremely attractive. In particular, an infrared diverging ghost dress- 



ing function is a necessary condition for the Kugo-Ojima construction 20 



Another is a non-perturbatively well-defined BRST symmetry [6(. Such a 



symmetry has been obtained [21], but is based on an average over Gribov 
copies. Thus subtle cancellation between Gribov copies is required to bring 
this construction into contact with the one presented here. 

If all of this would be correct in the end, this would be almost too good to 
be true: This would permit to use the gauge freedom to choose a gauge with 
useful properties. E. g., the explicit construction of the Hilbert space of QCD 
could be addressed by the use of the scaling (infinite B) Landau- B gauge 



and the Kugo-Ojima construction 20], which could provide one explanation 



of the confinement mechanism. Of course, also in the other Landau- B gauges 
gluon [sl, [iil and quark 23] confinement is obtained, though the explanation 



will necessarily be different. 
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